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I. INTRODUCTION 

It Is anticipated that Space Shuttle launches from Vandenberg Air 
Force Base (VAFBl. California, will be constrained to minimize the prob- 
ability of excee<. 9 design maximum wind loads. The reduction of wind 
loads Is normally achieved by biasing the ascent trajectory with respect to 
a nominal wind profile. The wind profile to be used and the rationale for 
its selection Is the subject of continuing studies. Ideally, maximum re- 
duction of wind loads would be achieved if wind biasing is with respect 
to the wind profile "seen" by the vehicle; such a profile is not observable 
prior to launch because of the variability of wind with respect to time 
and the inability of available wind measurement systems to measure winds 
in the space-time coordinate system of the vehicle. The time required for 
implementation of modifications to the wind bias program based on updated 
wind profile data also dictates that the wind profile measurement must be 
obtained prior to launch. 

The specification of an optimum time for the wind profile measure- 
ment depends in part on quantitative information on the temporal varii 
ability of winds aloft and knowledge of the load reduction achieved by 
utilizing wind profiles obtained at various time intervals prior to 
launch. Establishment of the optimum time interval is the objective of 
proposed simulation studies. A simulation study could Involve the use of 
a series of wind profiles measured at closely spaced time intervals. 

The last profile of the series is chosen to represent the profile seen by the 
vehicle and each proceeding profile is used for establishing a series of 
wind bias trajectories. The loads are calculated for each biased trajectory, 
and the amount of load reduction achieved is tabulated as a function of 
time interval prior to launch. The results are also compared with the load 
reduction that would have been achieved by simply using the monthly mean 
wind profile as the basis for the biased trajectory. An optimum time inter- 
val is chosen based on a statistical analysis of a large number of simula- 
tions. The ultimate implementation of the results of such a study is the 
establishment of a wind monitoring scheme to support space shuttle launches. 
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Quantitative Information on wind change with respect to time Is 
the basis for the simulation studies described above. This Is the second 
report prepared by SAI on the analysis of wind change at Space Snuttle 
launcn sites. This report is concerned with the Vandenberg AFB, California, 
launch site. The previous report llj was concerned with the Cape Kennedy, 
Florida, launch site. 

This report describes wind change at Vandenberg AFB for time 
intervals from 12 to 72 hours at altitudes from 0 to 27 km calculated 
from a ten-year (1965-74) sample of twice daily Rawinsonde profiles. Wind 
change with respect to time is presented in terms of statistical suninaries 
of wind component change and the modulus of vector wind change; the para- 
meters of idealized theoretical probability distribution functions 
representing the wind change variables are also presented. The validity 
of the hypothetical distributions is established by comparing them with 
the observed distributions. These distribution functions can be utilized 
to obtain statistical predictions of wind change with respect to time. 

Succeeding sections of this report consist of a brief statement 
of technical background (Section II), an analysis of wind change statistics 
(Section III), and conclusions (Section IV); the calculated statistics of 
wind change with respect to time by reference month at 1 km altitude 
intervals from 0 to 27 km are listed in the Appendix. 
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II. TECHNICAL BACKGROUND 

A. INTRODUCTION 

The large sample of wind profiles obtained at VAFB is suitable :\ r 
calculation of an equally large sample of wind change data. In ordr.' to 
I'eadlly abstract information on wind change from these data, it is 
necessary to perform a second series of calculations which provide 
statistical summaries of wind change. The choice of statistical parameters 
for description of wind change is based in part on the neeo to specify the 
parameters of theoretical distributions of wind change. These theoretical 
distribution functions are described in detail by Smith i2]. The basic 
distribution of the four va**iab1es consisting of tno zonal and meridional 
components of the vector wind at an initial time and after an elapsed time. 
At, is quadravariate normal. The conditional distnbiition c' the wind 
components at a specified future time, given the wind components at an 
initial time, is bivariate normal. The modulus of the wind change vector 
is Rayleigh and the distribution of either the zonal or meridional wind 
component change Is univariate noimial. A significant portion of the 
analytical discussion in Section III of this report is the presentation 
of observed distributions of wind change and comparison with the theoretical 
distributions of wind change variables. Succeeding paragraphs of this 
section are concerned with a description of the wind profMe data and 
the definition of statistical paranvters of wind change used in the 
various theoretical distribution functions. 

8. DATA 
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Wind change statistics for periods from to hours are calculated 
from the serially complete lO-:? kml Rawinsonde data at I km altitude 
intervals obtain»»u twice daily during the period 1965-T4 at Vandorberg 


AFB. California. The basic winds aloft data are recorded it. terms of 
wind direction, D and magnitude, W. The wind vector is expressed in the 
standard meteorological coordinate system m which the direction »rom which 
the wind is blowing is measured in degrees clockwise f**oi’i true nortn The 
zonal component, u, of the wind vector is positive for a west [wobt to 
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east) wind (9»270°) and negative for an east (east to west) wind (0=90°); 
the meridional component, v, is positive for a south (south to north) wind 
(0=180°) and negative for a north (north to south) wind (0=0°); u and v 
are obtained from 0 and W according to: 


u = -Wsin0, 360° 

(1) 

V = -W cos 0, 

(2) 

The relation between 0 defined above and the angle defined in 
standard mathematical polar form is: 


0 = 270 

(3) 

C. DEFINITIONS 


The subscript 0 is used to denote the initial value of a variable, 
and the suuscript 1 denotes the variable after an elapsed time. At. Thus: 

iu = 111 - 

(4) 
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Where Au and Av are the components of the wind change for a specified 
At, the modulus, R, of the wind change with respect to time is given by: 

R - ^ ( 6 ) 

The statistical means are denoted by an overbar, the standard de- 
viations and the correlation coefficients are denoted by and R(X,Y), 
respectively, vvith X and Y replaced with the notation appropriate to the 
variable of interest. 

0. STATISTICS 

The wind vector measurements at an ’niti-1 time and after an elapsed 
time are treated in this investigation as a sample from a quadravariate 
normal distribution defined by the fourteen statistics listed below: 
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The fourteen statistics of f’o k^jatirjvariate nornal dJStrihution v'' 
vector wind difference witn t-espect to ti"e •onsist of the fito h'vai’iate 


nomal statistics of vector wind at an iritia’ tiive ,j.. \ o o and 
'^0^' statistics 'tnolvinc comoonent ..1i fferentes wSion 


can be calculated f^'cm the duadrava'-'ate statistics listed above accordno 
to the following eoiiations: 
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is obtained from Equation 9 


V '('o- ’i' - V 


R(Vj. AV) . -L 


Av 


is obtained from Equation 10 
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in. ANALYSIS 


A. INTRODUCTION 

The statistics presented in the appendix of this report can be useful 
in the establishnwnt of a basis for certain aspects of Space Shuttle Launch 
planning for Vandenberg AFB. A prelaunch wind monitoring program may be 
required to provide data for assessment or modification of the Space Shuttle 
wind bias program. The development and utilization of the wind monitoring 
program will require knowledge oi^ the magnitude of vector wind change with 
respect to time. The analysis presented In this section establishes a 
theoretical basis for estimation of wind change. This is accomplished by 
comparison of theoretical probability distributions, whtch contain wind 
change sample statistics as parameters (from the appendix of this report ) . 
to observed probability distributions of wind change. Wind change with 
respect to time is analyzed herein in teir.is of wind component change, un- 
conditional and conditional joint distribution of wind component change 
and the modulus of vector wind change. 

B. WIND COMPONCNT CHANGL WITH Rf.SI'LC’ VO HMl 

The theoretical probability distribution of wind component change 
with respect to time is univariate normal with zero moan and standard 
deviation given by Fguations and 10; the assumption of zero means of 
component differences is verified by the sample statistics given tn the 
appendix. The theoretical nonuxl distribution of component differences can 
be derived by using either the standard deviations of component difft'rences 
givett in the appeodix or an estimate which c.tn be obtained from the staiul.xrd 
deviation of tlio components if it is .tssumed that: 
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Equations 9 and 10 reduce to 


=V^ Oy V^l - R(uj. Uq) 


°Av ‘V^ ^v Vl - R(vi» ''o^ 


The wind component autocorrelation functions, R(uj, Uq) and R{Vj, Vq) can 
be represented by a negative exponential function of time increment, 

T j i « e . , 

R{uj, Uq) • EXP (-bt) {U 

R(v^, Vq) • EXP (-CT) (IS 

where b and c are computed according to 
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Examples of the decay of the autocorrelation function at 12 km during 
November, December and January at VAFB are illustrated in Figure 1; the 
lines in the figure represent the decay rate predicted by Equations 18 and 
19. 

Substitution of Equations 18 and 19 into 16 and 17, respectively, 
yields a simple expression for and in terns of and o^, 
respectively. 
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Figure 1. Zonal and Meridional Wind Component Autocorrelation 
at 12 km at VAFB (1965-74) 
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Equations 22 and 23 Indicate that and are asymptotic to 
y[z Oy and a^ for large values of t. Therefore, estimates of the 
extreme value of and are obtained by setting t equal to “ In 
equations 22 and 23. 

The calculated values of b and c for VAFB during November, December 
and January listed In Table 1 are also plotted In Figures 2 through 4. 
Disregarding minor variations at small altitude Intervals, the decay of the 
zonal component autocorrelation function with respect to time decreases with 
altitude during November, December and January. The decay of the meridional 
component autocorrelation function generally exceeds the decay for the 
zonal component during November, at altitudes above 16 km during December 
and above 1 km during July. The most rapid decay of the autocorrelation 
function occurs for the meridional wind component at altitudes between 8 
and 16 km during January. The calculated and observed values of 
and o. (t) at 1, 6, 12, 18 and 24 km during November, December and January 
are listed in Tables 2 through 4. The estimated extreme values of and 
o.„, a and VFo„. respectively), are listed at the bottom of each 

4jV U V 

column of calculated values. The comparisons In Tables 1 through 3 indicate 
that 0 ^^ and 0 ,^^ can be accurately estimated by application of Equations 22 
and 23, respectively. General application of this estimation technique at 
other locc;,ions utilizing published statistics of wind component standard 
deviations (as in [3], for example) would require a more adequate knowledge 
of the form of the autocorrelation function than Is presently available. 


The theoretical distribution of wind component differences has been 
derived from sample estimates of and and lu and Av (given in the 
appendix) fc the intervals of 12, 24, 36 and 48 hours during November, 
December .nd January at 12 km over VAFB; the theoretical normal distributions 
are plotted as straight lines in Figures 5 through 10; the plotted symbols 
represent the observed distributions of Au and Av. It is indicated that the 
observed distribution of component changes is accurately represented by the 
theoretical normal distribution for a large range of probabilities; the 
deviation of the observed distribution from the theoretical distribution at 
the extreme probabilities Is attributed to the small sample of observations 
ai.a errors in the Rawinsonde data. 
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Table 4. Calculated (Eqs. 22, 23] and Observed a.„ and a. During 
January (1965-74) at VAFB 
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Figure 6. Theoretical (straight lines) and Observed (Plotted Points) 
Cumulative Probability Distribution of Meridional Wind 
Component Change, ^v, with respect to time Increment, t. 
During November at 12 km at VAFB (1965-74) 


i 


{ 



f 


j 


} ’ 


I HlfilNAL PA(JK 18 

"f'' guAU'ry 


N 













■:v'';''Tn"^rrrTi-^T:r.TZT: 





7 ■- • 


< 

) 



Figure 9. Theoretical (straight lines) and Observed Plotted Points) f ■ 

Cumulative Probability Distribution of Zonal Wind Component 1 ! » 

Change, au, with Respect to Time Increment, t. During January if ' 

at 12 km at VAFB (1965-74) I' 
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C. JOINT DISTRIBUTION OF WIND COMPONENT CHANGES WITH RESPECT pn TTMF 

The joint distribution of zonal and meridional wind component change 
with respect to time (Au and Av) can be approximated by a bivariate normal 
distribution. A useful property of such a distribtion is that an ellipse 
can be calculated which contains the end points of a specified percent of 
vectors having components Au and Av. A detailed description of the de- 
rivation of probability ellipses and plotting methodology is aiv»en by 
Smith [2]. The five parameters of the bivariate normal distribution of 
Au and Av, calculated for each monthly reference period at VAFB at 1 km 
altitude internals from 0 to 27 km are listed in the apoendix. 

The degree of approximation of the bivariate normal distribution to 
the observed distribution can be evaluated by comparison of the observed 
percentage of vectors which are contained within the ellipse to that pre- 
dicted by the ellipse at a specified probability level. For example, for 
a sample of 1,000 vectors, 950 of the vectors should terminate within the 
95 percent (theoretical P = .95) ellipse calculated from the bivariate 
statistics of the 1,000 vectors; however, a plot of the 1,000 vectors 
could indicate that only 945 vectors (observed p=.945) terminate within the 
95 percent ellipse. For illustration on a linear graph comparison of the 
theoretical to the observed P is given in terms of the parameter X given 
by 

Ag = yfz V- In (1"P) (24) 

A comparison of theoretical and observed values of X^ for November, 
December and January at 12 km for time intervals of 12, 24, 36 and 48 
hours is illustrated in Figures 11 through 13. Perfect agreement between 
theoretical and observed X^ is represented by a line drawn from the origin 
with a slope, B, equal to 1. The calculated least squares slopes are 
given in the figure legend. The plots indicate a tendency for the 
theoretical X^ to exceed the observed \„ for large values of \^. The 
interpretation of this result is that for extreme probabilities the theon - 
ical distributions predict fewer wind change vectors terminating outside 
the ellipse than is observed. These results may have to be taken into 
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V 

consideration it engineering application of theoretical wind change 
statistics beyond the 95 percent level is required. 

The 95 percent probability ellipses for the joint distribution of 
wind component changes with respect to time at 6, 12, 18 and 24 km during 
November, December and January are illustrated in Figure 14; the relatively 
small changes with respect to time at 18 and 24 km, the large changes at 
12 km and the similarities between the three months are clearly il lusti'ate.l. 
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Observed \ as a Function of Theoretical X for a Bivariate 
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Respect to Time at 12 km During December 
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D. 


MODULUS OF VECTOR WIND CHANGE WITH RESPECT TO TIME 


If wind change with respect to time has a distribution which is 
bivariate normal, the modulus R, of the wind change vector (defined by 
Equation 6) has a Rayleigh distribution. Since the Rayleigh distribution 
cannot be integrated in closed form, numerical integration is required to 
obtain the cumulative probability distribution. Derivation of the Rayleigh 
distribution, given the five bivariate normal distribution statistics, 
requires summation involving products of the modified Bessel function of 
the first kind. Smith [2l summarizes the basic equations for the Rayleigh 
distribution derived by Wier [4] and extended by Yadavalli [5] to include 
the condition for correlated variables. The Rayleigh distribution reduces 
to the integrable classical forr. if it is assumed that the components of 
the vector wind change are independent and that they have zero means and 
equal standard deviations; the classical Rayleigh probability density 
function is 

f(R) » \ EXP (-R^/2o^) R = 0 (2b) 

Integration of Equation 25 from zero to a specified value of R yields 
the cumulative probability that R ^ R* where, 

Pr |R V R*| = 1 -EXP (-R^/2o^) R = 0 (26) 

where c- ^ = e. 

.iu .\v 

Since the standard deviation of the component difference can be 
expressed as a function of the standard deviation of the components 
(Equations 22 and 23) it follows that 


Pr |r < R*) = 1 - EXP 

where \ and k correspond to either o 
k u 

22 and 23. 




1 1-EXP (-k-. )1 


(27) 


and b or and C given in Equations 
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An expression for R given a particular probability, Pr [R < R*i , is 
obtained by solution of Equation 27 to obtain 


R = V 2 ^e^ Vl - EXP (-kr) 
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where X^ is derived from Equation 24 denoting Pr (R v R*j by P 

The choice of = 0 ^ and k * c (from Equation 23) at 12 km during 
November, December and January yields the most accurate approximation of 
the cumulative Rayleigh distribution obtained by numerical integration of 
Equation 28 in Reference 1. Comparisons of the 99, 95, and 50 percentile 
modulus of the wind change vector with respect to time based on the Rayleigh 
(Equation 28, Reference 1) and the classical Rayleigh (Equation 27) are 
illustrated in Figures 15, 16 and 17; the rather good agreement between 
the distributions for time intervals from 12 to 72 hours is attributable 
to the accuracy of the simplifying assumptions described above. There 
is a slight tendency, especially for time intervals < 36 hours for the 
classical Rayleigh to be larger than the Rayleiqh. 

The remaining question is: How well do these theoretical distribut iv>ns 

compare with observed distributions? Comparisons of observed and theoretical 
values of R for time intervals of 12, 24, 36 and 48 hours at 12 km during 
November, December and January at VAFB are given in Tables 5 through 7; 
column II of the tables contains R calculated according to the classical 
Rayleigh distribution with 0 equal to the monthly value of at 12 km and k 
equal to the decay constant, c. in the monthly exponential least squares 
fit to the V component autocorrelation function (Equation 23); Column I 
was obtained b\ numerical integration of the Rayleigh distribution. It 
is indicated t!;at the observed cumulative distribution agrees fairly well 
with the theoretical distribution for probabilities less than .96 to .97. 

For large probabilities, .nere is a consistent tendency for the theoretical 
distribution tvi underestimate Che observed di stribution. This tendency is 
attributable to the small samole of data available at the extreme probabil- 
ities and errors in the Rawinsonde data. 
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Table 5. Theoretical and Observed Modulus, R, of Vector Wind 
Change with Respect to Time at 12 km during November 
(1965-74) at VAFB 
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f. CONDITIONAL VtCTOR WIND LLLIPSES j' 

! V 

I * 

Prior knowledge that environmental constraints necessary to assure 1 

the success of a space vehicle launch will be satisfied implies that there 
is a capability for prediction of environmental parameters; the predictKri 
can be based on knowledge of conditions prior to launch. With regard to 
winds aloft, prior conditions are typically based on Rawinsonde or Jimspheie 
wind profiles. A typical question that could be posed before launch is: I 

Given a measurement of the wind vector \2 hours prior to launch at km. 
will the wind vector at launch time be within 95 percent reference month 
wind ellipse? A question ot this typo can be answered if the disti ibutien 
of vector wind components at an initial time, T^^. and at a future time. 

Tj. can be approximated by a quadravariate normal distribution. Gi\en Nie 
components of the vector at 1\^, the conditional distribution of the vect.'; j 

wind at T. IS bivariate normal. Smith [!i vlescribes the derviation of 
the conditional bivariate nonnal distribution and documents the computer 
program used in this investigation for calculation of these distributions. 

Figures lb thru .'0 illustrate the d*; percent conditional bivariate 'ormal 
distributions at IJ km that have been calculated for time increments of j 

L', L'4. 3b. 48, bO and 7;? hours for the months of November, December and 
January. Five vectors were selected as given initial conditions for 
calculations of the conditional ellipses. Ihe components of the vectors 
aro defined below. 

1. Voiit'ilv componont moans given by Falls ;.li. 

M.iximuiii conal wind and the corresponding aorulional wind fro'i 
; !ie monthly d;, percent vector wind ellipse. 

'. Mimnuim .:onal wind and the corresponding meridional wind from 
t '10 monthh dp percent vector wind ellipse. 

4. ^'a\^mu!ll meridional wind and the corresponding zonal wind from 
the monthl.v bb po.-cent vector wind ellipse. 

Ml niimin meridional wind and the correspending zonal wind i rom 
the monthly dp percent vector wing ellipse. 

The viiven vectors are specified in the inset of Figures 18 through .)0 j 

lP''1ar term, at kml and in Tables 8 througn 10 tcomponont form, at 1. ' 

18 an.l .'4 kmV 
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The conditional ellipses illustrated at the center of Figures 18 
through 20 show that if the observed wind vector has components equivalent 
to the monthly mean components (condition 1) then 95 percent of the wind 
vectors after elapsed times as large as 72 hours will fall within the monthly 
95 percent ellipse. Therefore, satisfaction of a launch constraint which 
states that the wind vector must be included within the 95 percent monthly 
ellipse would be assured for periods as long as 72 hours following an obser- 
vation of a wind vector having components which correspond to the monthly 
means. The conditional ellipses based on selection of given wind vectors 
that terminate on the monthly 95 percent ellipse (conditions 2 through 5) 
have a significant proportion of their area lying outside the monthly 95 
percent ellipse; as the time increment increases this proportion decreases 
but remains significant for a time increment as large as 72 hours. This 
implies that a significant proportion of wind vectors will not satisfy a 
launch constraint based on the 95 percent wind ellipse for periods as long 
as 72 hours (or longer if these calculations are extended) following an 
observation of a wind vector which terminates on the 95 percent ellipse. 

The wind direction characteristics of a wind ellipse can be described 
in terms of the angles associated with wind vectors constructed between the 
origin and the center of the ellipse (at the component means) and between 
the origin and the two tangent positions to the ellipse. The three vectors 
constructed in this manner and the angles 0£ and A6 are illustrated 

in Figure 21, the range of wind angles, 9j^, is 9^ to 9g. The angles 9^, 

AO calculated from five 95 percent conditional ellipses for November, 
December and January at 1, 6, 12, 18 and 24 km are listed in Tables 8 
through 10. 
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Table 


8. Wind Direction (Degrees) Characteristics of 95 Percent 

Conditional Wind Ellipses at VAFB during November (1965-74) 
for an Elapsed Time (t) of 12 hours 
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Table 9. Wind Direction (Decrees) Characteristics of 95 Percent 

Conditional Wind Ellipses at VAFB during December (1965-74) 
for an Elapsed Time (r) of 12 hours 
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Table 


. Wind Direction (Deorees) iharactenstics of 05 PerctMU 
Conditional Wind ellipses at VAFB during January (1965-74) 
an riapsod Time (; ) of doufs 
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YEAR TO YEAR VARIABILITY 

Artificial wind profiles that are constructed using refe»'ence wnth 
statistics, such as those based on an eight-year sample given by Falls (3), 
have been used for development of wind bias programs for launch vehicles. It 
is of interest to study how realistic these reference month statistics are 
when compared with data from individual months. 

A comparison of individual monthly means to reference month me^ns given 
by Falls [3] is illustrated in Figures 32 and 23 for the zonal and meridional 
component, respectively. The reference month means are illustrated as a 
repeating function. The distribution of deviations of zonal and meridional 
wind components from the reference month means illustrated in Figures 23 and 
.13 can be treated as havina a normal distribution. The distribution of the 
deviations during the period 1965-74 are illustrated in Figure 24; the 
straight lines represent normal distributions and the plotted symbols re- 
present the observed distribution. It is indicated that 10 percent of the 
monthly means at 12 km have an absolute difference of more then S.5 m/sec 
from the reference month means. 

Another indicator of the year to year variability is the variability 
of the bivariate normal statistics of tne vector wind. These statistics 
at 12 km for any months during the period 1965-74 are listed in Table 11. To 
illustrate the variability, tne 95 percent probability ellipse for the wind 
at 12 km during Novemoer 1965 and 1969 are shown in Figure 25. 

"he variability of the wind statistics from year to y< '.r described 
above supports the conclusion that the reference ironth statistics based on 
a relatively long data sample do not necessarily provide an accurate estimate 
.'f the conditions during a particular month. If wind biasing is to be 
required for reduction of wind loading on the Space Shuttle, the utilization 
of an artificial wind profile composed of reference month sampling period 
winds ’ay r.ot vield the desired results. It is ’'ecommended that pre-launch 
wind bias models be based on sots of wind profiles that moresent the range 
of year to year variability indicated nereiti. In .in o;vrational sense, the 
greatest reduction of in-flioht loads would be achieved b\ wind biasing with 
respect to a wind profile measurement taken prior to launch. 
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Bivariate Normal Statistics of the Vector 
Wind at 12 km Over VAFB for Individual Months 
During the Period 1965-74 
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Bivariate Normal Statistics of the Vector 
Wind at 12 km Over VAFB for Individual Months 
During the Period 1965-74 (Cont'd) 
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Figure 25. 95 Percent Vector Wind Ellipse at 12 km During 

November 1965 and 1969 at VAFB 
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G. COMPARISON OF WIND CHANGE WITH RESPECT TO TIME AT KSC ANO VAIH 

Wind change at KSC has been studied in detail in an earliei ixn'ort il; 
In anticipation of questions that may arise concerniiiq comparisons of Kf.C 
and VAFB, a set of 95 percent vector wind change ellipses at 12 km ha\e 
been prepared for the two locations. The ellipses illustrated in I loure -’v 
are for the months of January, April, July and November. The Ijrqei' wind 
change at VAFB is clearly shown during all months with the exception of Juli 
when the wind change is small at both locations. 



JULY 




APRIL 




JANUARY 




*THE SIX ELLIPSES FOR EACH MONTH ARE FOR THE INCREMENTS AT 12 KM HR INTERVALS 
FROM 12 TO 72 HOURS: THE AREA OF THE ELLIPSES INCREASES WITH INCREASING 
TIME INTERVAL 

Figure 26. Comparison of 95 Percent Wind Change Ellipses at. 12 km for 
VAFB (1965-74) and KSC (1956-70) 
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IV. CONCLUSIONS 

The analysis presented in the preceding section for selected months 
and altitudes illustrates how various theoretical distribution functions 
can be used for c. 'ulation of wind change with respect to time at 
Vandenberg AFB, California. The calculations can be made by utilization of 
the statistics given in the appendix for any reference month at 1 km 
altitude increments from 0 to 27 km. 

The basic underlying assumption for the calculation of the distribu- 
tions is that the joint distribution of the four variables represented by 
the components of the wind vector at any initial time and after a specified 
elapsed time is quadravariate normal. If the wind vector is specified at 
an initial time, then the conditional joint distribution of the wind com- 
ponents at a future time is bivariate normal. Since each of the variables 
of the quadravariate normal distribution is normal and the difference of 
two normal distributions is normal, it follows that wind component change is 
also normal and the joint distribution of zonal and meridional wind change 
is bivariate normal. The modulus of bivariate normally distributed" variables 
has a Rayleigh distribution, therefore, the modulus of vector wind change 
with respect to time is Rayleigh. 

Sample distributions based on reference month Rawinsonde data obtained 
during 1965-74 agree reasonably well with the aforementioned theoretical 
distributions. 

The standard deviation of wind component change with respect to time 
is the only statistic required for determination of the theoretical pro- 
bability distribution (normal with zero mean) of wind component change. 
has been shown that over a large range of altitudes that this statistic 
can be estimated from wind component standard deviation and the decay 
constant of the component theoretical autocorrelation function (Table 1). 

The assumption of exponential decay of the autocor*'elation function is 
reasonably accurate in most instances to time increments as large as 60 
hours during November, December and January. 
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The observed modulus of vector wind change with respect to time Is 
systemmatlcally larger than the predicted modulus (Section III.C) for 
extreme probabilities. This may be attributable to inadequacy of the 
theory or Inaccuracies of the data which affect the observed distribution 
at the extreme probabilities. If the theoretical distribution at extreme 
probabilities Is to be used In engineering applications. It will be 
necessary to explain these systemmatic differences. 
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I APPENDIX 

This appendix contains two sets of reference month quadravarlate and 
conditional bivariate normal statistics of variables X, Y, XP and YP, at 

1 

j 1 km Intervals from 0 to 27 km. The statistics were calculated from IG 

^ years (1965-74) of twice dally VAFB serially complete Rawinsonde data. The 

notation for the variable given in Section II of this report differc from 
! the notation established for the computer output given herein; the 

notations are compared In Table A-1. 


I I 


) . 


TABLE A-1. NOTATION OF VARIABLES 


Computation Set 

A 

B 1 

Variable 

.. - . . 

Text 

(Sect. II) 

Computer 

Output 

Text 

(Sect. II) 

Computer 

Output 

X 

^0 

u(at T) 

^0 

u(at T) 

Y 

''o 

v(at T) 

''o 

v(at T) 1 

1 

' XP 

1 
1 

“l 

u(at T+DT) 

Ui - Uq 

= Au 

u(at T+DT) I 
-u(at T) 

\ 

YP 

1 

1 

L 

''1 

v(at T+DT) 

0 
> 

1 > 

> II 

v(at T+DT) * 
-v(at T) j 

I 

1 


Table A-1 shows that the quadravarlate statistics of computation set 'A" 
are for wind components at an initial time and after a specified time 
increment; the statistics for set "B" are for wind components at an initial 
time and wind component change after a specified time increment. The re- 
ference month quadravarlate normal statistics at a particular altitude for 
six time increments (12, 24, 36, 48, 60 and 72 hours) are listed in the 
lower left of each page of computer listing; the six sets of conditional 
bivariate normal statistics corresponding to the six time increments are 
listed in the lower right. The data were conditioned monthly means 
given by Falls (3l. The derivation of the conditional bivariate normal 
statistics for any other given vector involves recall elation of the 
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conditional means according to equations A-l and A-2; the standard de- 
viations and correlation coefficients do not have to be recalculated 
because they are independent of the given wind vector. 

[(R(x,xp) - R(x,yp) R(xp,yp)) (xp* - 5^) 

+ (R(x,yp) - R(x,xp) R(xp,yp)) (yp* - yp) (o/c)} 

X |xp*=x + j (A-n 

1 - [R(xp,yp)j^ 


t(R(y,xp) - R(y,yp) R(xp,yp)) (xp* - xp) (« 

_ _ + (R(y.yp) " R(y,xp) R(xp,yp)) (yp* - fp") (a ,, ... 

y |yp.=y . , lA-., 

1 - fR(xp,yp)]^ 

where, x and y are the mean components of the conditional distribution, 
xp* and yp* are the components of the given vector and 

a ,0 ,0 and o are equivalent to S.D.x, S.D.y, S.O.xp and ! 

X y xp yp 

S.D.yp, respectively given in the computer listings. 
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